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I. INTRODUCTION 

Among the so-called 'high precision' nucleon-nucleon (NN) potentials the Argonne V18 potential [l[ is the only one 
given in a basis of spin-isospin operators, O n , which are multiplied by scalar functions, V n (r), of the relative distance 
r. These features are advantageous for some applications, however, there are classes of problems where a momentum- 
space treatment is preferred. These applications include rclativistic fcw-nuclcon scattering and electroweak probes of 
few-nucleon systems. While realistic momentum-space nucleon-nucleon potentials are available, they are either given 
in terms of partial- wave expansions, e.g. the CD-Bonn interaction and the Nijmegen interaction Q, or they are 
limited to low energies, like the interactions derived in chiral effective field theory [H, 

The experience in standard thrce-nuclcon calculations that is based on a partial-wave projected momentum-space 
basis @ shows that while this standard treatment is quite successful at lower energies, the numerical realization of 
the scattering equations becomes more tedious with increasing energy. For a system of three bosons interacting via 
scalar forces it already has been demonstrated that it is relatively easy to directlycalculate three-body scattering 
observablcs in a relativistic Faddeev scheme without using partial wave projections [3, Thus, it is natural to strive 
for solving the three-nucleon Faddeev equations in a similar fashion. 

The first step in this direction is solving the two-nucleon problem based on an operator expansion instead of a 
partial-wave projected momentum-space basis. One of the advantages of using operator expansions over partial waves 
is that the transition matrix elements are reasonably smooth. At higher energies, an accurate representation of these 
smooth, but forward-peaked, functions requires a large number of partial waves. A converged partial-wave expansion 
of the transition matrix elements implies that all of the high-frequency oscillations from the large-/ partial waves must 
largely cancel. In addition, to achieve these cancellations the computation of the transition matrix elements requires 
the accurate computation of oscillatory integrals. At low-energies partial waves are preferred because they replace the 
two-variable Lippmann-Schwinger equation by a finite set of uncoupled one-variable equations. For sufficiently high 
energies the simplicity of a direct vector treatment of the Lippmann-Schwinger equation has advantages, even though 
the integral equation has an additional integration variable. In this paper we discuss the direct vector solution of the 
Lippmann-Schwinger equation using a recent momentum-space treatment of the Argonne V18 potential There 
the potential was only tested in the calculation of the deuteron bound state, which is not sensitive to the charge 
symmetry breaking parts of the AV18 potential. Solving for the scattering observables provides a complete test of 
the momentum-space representation of the AV18 potential presented in Refs. |9l. Il0|. 

There have been several approaches formulating nucleon-nucleon (NN) scattering without employing a partial- wave 
decomposition. A helicity formulation related to the total NN spin was proposed in [ll|. The spectator equation for 
relativistic NN scattering has been successfully solved in [l2j], also using a helicity formulation. 

A three-dimensional formulation based on an operator expansion was proposed and carried out for a chiral next-to- 
ncxt-lcading order NN force and a standard one-boson-exchange potential in [l3[. A basic foundation for the latter 
rests on the fact that the most general form of the NN interaction can be represented as a linear combination of 
six linearly-independent spin-momentum operators with scalar coefficient functions. This representation determines 
the spin-structure of the NN bound and scattering states. The Wolfenstein decomposition of the NN scattering 
amplitude into five linear operators is dictated by rotational, space reflection, spin-exchange symmetry, and time- 
reversal invariance (T3 |. A sixth independent operator with these symmetries does not exist on shell. Ref. [l3j uses 
six linearly independent operators that satisfy all of the symmetry requirements, but these operators become linearly 
dependent on shell so the scattering amplitude was computed from the off-shell result by the required continuity of 
the transition matrix elements. 

The vector treatment of the spins requires the computation of an analytic expression for each spin-basis element 
at each quadrature point. A symbolic reduction technique, developed in (Toj . automates the computation of these 
expressions. Our calculations use five independent spin-operators to expand the on-shell transition matrix element and 
an additional sixth operator to expand the off-shell transition matrix elements. Based on calculations using different 
choices of the sixth operator, we find that using a sixth operator that changes sign under time reversal to expand 
the off-shell transition matrix elements leads to a numerically more stable discretization of the Lippmann-Schwinger 
equation when compared to calculations using a sixth operator that has the symmetries of the potential. An important 
observation is that while a time-odd operator appears in the expansion of the potential, the potential itself remains 
invariant with respect to time reversal because the expansion coefficients also contain time-odd components. The 
expansion coefficient of the time-odd operator that we consider vanishes on-shell. 

In Sec. II we discuss the operator basis that we use to expand the NN potential and derive in Sec. Ill the resulting 
Lippmann-Schwinger equation and show how we explicitly obtain the on-shell scattering amplitude. We then discuss 
our numerical procedure and the extraction of the Wolfenstein amplitudes in Sec. IV. In Sec. V we discuss the selected 
NN observablcs at different energies, and conclude in Sec. VI. 
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II. OPERATOR EXPANSIONS 



Momentum-space scattering calculations are performed using the Lippmann-Schwingcr integral equation. It is a 
singular integral equation with a compact kernel. This means that the kernel can be uniformly approximated by a 
finite dimensional matrix and that the equations can be solved to any desired accuracy by solving a finite system of 
linear equations. 

The Lippmann-Schwingcr equation has the well-known form 

T(z)^V + V(z-H )- 1 T(z). (1) 

We represent the operators in this equation by matrix elements in the total momentum P and single-particle momen- 
tum, p = 7}(p P — Pn)j boosted to the system rest frame with a Galilean boost. Matrix element of the interaction and 
transition operator have the form 

<p',p'|F|p,p> = <5(p'-P)<p'||y|| P ) (2) 

and 

(P\p'|T(*)|P,p) = S(P' - P)(p'||T(z)||p), (3) 

where the reduced matrix elements (p'||^||p) and (p'||T(z)||p) are matrices in spin-isospin space. 

Any operator V on the two-nucleon Hilbert space can be expanded as a linear combination of the spin-basis 
operators, £ M „ := (g> a 2lJ , where a ifJ/ = (I, a), 



V = J2 (4) 



and 



V" = -Tr(^V). (5) 

The operator E^j, has sixteen components, however the number of independent operator types that can appear in 
a transition operator or potential is constrained by rotational invariance, time-reversal invariancc, space-reflection 
symmetry, and spin-exchange symmetry. Spin-exchange symmetry means that the potential commutes with the 
square of the total spin, which is a symmetry of nuclcon-nucleon interactions. These symmetries can be utilized to 
reduce number of coupled Lippmann-Schwingcr equations. 

The three Pauli spin matrices in E M „ for each particle can be replaced by the rotationally invariant operators, 
Vj • Oj, for any independent set of vector operators, {Vi, V2, V3}. The traces of the product of these operators with 
the potential have the same time-reversal, space-reflection, and spin-exchange symmetry as the operators. If the spin 
operators are chosen to be orthonormal with respect to the trace norm, the only operators that contribute to the 
potential expansion are the operators that satisfy all these symmetries, except time-reversal invariance, because it 
is possible to have a scalar coefficient function that changes sign on time reversal. Consideration of the symmetries 
of a complete set of spin operators implies that the most general interaction consistent with these symmetries can 
be represented by an expansion in terms of five operators on shell and six off shell. When the spin operators are 
not orthonormal with respect to the trace, a matrix inversion is needed to compute the expansion of the potential. 
This inverse matrix may not have any simple symmetry with respect to time reversal. The same will be true for the 
resulting expansion coefficients. Nevertheless, the number of operators needed is unchanged and the symmetries of 
the full potential are preserved. 

In this paper we consider the following independent momentum operators used by Wolfenstein [l5| 

K 

Q 

N:= 



(p' 


-p) 


IP' 


-Pi 


(p' 


+p) 


IP' 


+Pl 


(p' 


x p) 


IP' 


x p| 



These vector operators can be classified by their transformation properties with respect to space reflection and time 
reversal. 
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A useful independent set of spin operators that arc sufficient to expand the scattering amplitude matrix M was 
given by Wolfenstein. In terms of the vectors of Eq. ^ the Wolfenstein operators are 



{/, K • Q) ® (<x 2 • Q), (a t ® 1 + 1 (8 <r 2 ) • N, (<n • N) ® (<r 2 • N), (<ti • K) g> (cr 2 ■ K) j. 



(7) 



The coefficients of the expansion of the scattering amplitude matrix in Wolfenstein spin-operators are the so-called 
Wolfenstein parameters, a, c, m, g, h. Knowledge of these parameters as a function of the initial and final momenta 
contains all of the information in the scattering amplitude matrix. The scattering-amplitude matrix is related to the 
on-shell transition matrix element by 

a/(p',pH-4^<p'II^)IIp> (8) 

where fi is the reduced mass of the two-nucleon system and the Wolfenstein parameterization of M(p', p) is 

M 



al + c(cti ® 7 + J <8> <r 2 ) • N + m(<n • N) <g> (er 2 • N)+ 
(g + h)(tn • Q) ® (<r 2 • Q) + (g - fyfa ■ K) ® (<x 2 • K). 



(9) 



The operators of Eq. (|9]) are invariant with respect to rotations, space reflection, spin exchange, and time reversal. 
The additional operator, 



Oi • K) ® (a 2 ■ Q) + (o-i • Q) ® (<t 2 • K) 



(10) 



is rotationally invariant, space reflection invariant and spin exchange invariant, but changes sign under time reversal. 
If this operator is paired with a coefficient function that also changes sign under time reversal then the product is 



invariant with respect to time reversal. A time-odd coefficient function must be odd in |p 



l\2 



and thus vanish 



on shell. This means that while this operator cannot appear in the scattering amplitude matrix, M, it can appear 
in the interaction or half-shell transition matrix and it may also appear in the off-shell unitarity constraint on the 
transition operator. 

The six independent operators of Eqs. (0 and (|10j) arc sufficient to expand any potential that is invariant with 
respect to rotations, space reflection, time reversal, and spin exchange. 

In order to solve the Lippmann-Schwinger equation the interaction is expanded in terms of these six independent 
spin operators. In our applications the operator (<Ti • N) ® (<r 2 • N) is replaced by the operator (eri • cr 2 ). This 
replacement is done for numerical reasons, but it is also useful because this operator appears in the spin-spin part of 
the interaction as well as in the tensor force. 

The transition matrix elements can be expanded using the same six operators. This results in a set of six coupled 
equations for the half-shell or off-shell transition matrix elements. Since the scattering amplitude matrix is related to 
the on-shell transition operator by a multiplicative scalar, it is useful to put the integral equation in a form where the 
on-shell points are included in the linear equation. To avoid ambiguities, the on-shell amplitude is expanded using 5 
independent spin components. This will be discussed in section IV. 

The expansion of both the interactions and transition matrix elements is constructed by taking traces with the six 
independent spin operators 



w 1 


= I 








w 2 


= <Tl • 








w 3 


= (ffl 


■K)« 


5 (<T2 


•K) 


w 4 


= K 


•Q)s 


5 (CT2 


•Q) 


w 5 


= i((Ti 


(81- 


-I® 


c 2 ) • N 


w 6 


= (<ri 


■K)« 


5 (ff2 


• Q) + (a-i 



Q)®(<r a -K). 



(11) 



It follows that 



V 



(12) 



where 



Ui := (A _1 ) y, Er(W J V) A 1 ? := Tr(W l W j ). 



(13) 



Both T(z) and V are matrices in the spin and isospin degrees of freedom and the and depend on p and p'. 
By taking traces, the Lippmann-Schwinger equation becomes a set of six coupled integral equations in two variables. 
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As discussed before, six linearly independent operators are required to expand any NN potential. This expansion 
can include the product of a time-odd, space reflection invariant spin-momentum operator with a space reflection 
invariant, time-odd coefficient function resulting in a term that conserves the symmetry properties of the NN-potential. 
Note that a space-reflection odd operator cannot appear in this expansion because it would require a pseudo-scalar 
coefficient which cannot be constructed from two vectors (Toj . 

The spin-spin operator, W 2 , and the Wolfenstein operator (avN)(a- 2 -N) are independent and have all of symmetries 
of the potential. They are related by 



R 2 K 2 Q 2 



(K 2 Q 2 - R 2 ) 2 (K 2 Q 2 - R 2 ) 2 



((<ri-K) (tT 2 -Q) + ( < 7 1 .Q)( < r 2 .K)) 



Q 2 R 2 K 2 Q 4 \ (a 1 -N)(cT 2 -N) 

(CTl ■ JVJ (<T 2 • KJ 



(K 2 Q 2 -R 2 ) 2 [K 2 Q 2 -R 2 ) 2 ) N 2 
K 2 R 2 K 4 Q 2 ^ 



(K 2 Q 2 - R 2 ) 2 (K 2 Q 2 - R 2 ) 2 



(o-i • Q) (<ra • Q) , (14) 



which involves W , W as well as the time-odd operator, W . In this expression R — K • Q = |p'|~ — |p| is the 
time-odd coefficient of the time-odd operator W 6 , which vanishes on-shell. 

It follows that the basis for the off-shell potential, W 1 — W 6 could be replaced by an equivalent basis where W 6 
is replaced by the Wolfenstein operator (<xi ■ N)(er 2 • N). While the second basis has all of the symmetries of the 
potential, the basis that includes the time-odd operator W e is preferable for numerical reasons as will be elaborated 
on later. 



III. LIPPMANN-SCHWINGER EQUATIONS 



The reduced transition matrix elements can be expressed as 

t(p',p,x;z) := (p'||T(z)||p) 



z:=f°+*0+, 



(15) 



where the variables are defined below and t(p' ,p, x\ z) is a spin-isospin valued function. Here po represents the on-shell 
momentum. The Lippmann-Schwingcr equation for t(p',p,x; z) is the two- variable integral equation 

t(p',P,x';z) = 



v(p',p,x') + / dx" / dp"p" 2 u(p',x',p",x")G (p",z)t(p",p,x";z) 



(16) 

where the coordinate system is chosen so that the initial momentum is in the 3-direction and the scattering plane is 
the 1-3 plane. The integration over the azimuthal angle only affects variables in the potential and can be carried out 
independently. The momenta and quantities in Eq. (|16p are parameterized in terms of the variables p,p',x',p",x" 
defined by 

p := (0,0, p) = p 



P := 



P P 

P"P 
P'p" 



v(p',p, x' 
v(p',p",y) 

u( P ',x> \ P " \ x " ;</>") 

G (p",z) 
u{p ,x ,p ,x ) 



(p'yjl -f 2 ,o, P y) 

(p" \fl - x" 2 cos(0") , p" \fl - x" 2 sin(0") , p"x" 

x' 
x" 

y = x'x" + ^(l-x ,2 )(l-x" 2 )cos((f)") 

V(p',p) 

V(p',p") 

v(p',p",x'x" + V '(1 - x' 2 )(l - x" 2 ) cos((t)")) 
pi - p" 2 + i0+ 

j ill I I I II II ±11 \ 

d<p u{p ,x,p ,x ,<p ). 
o 



(17) 
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In the following we omit the explicit limits of the integrals. One can see that Eq. (fT6|) is an integral equation in two 
variables, the integration over the azimuthal angle <f>" can be carried out independently as given in Eq. (|17[) . 
In order to reduce Eq. to an algebraic equation the singular integral is first treated by a subtraction 

t(p',p,x';z) = v(p',p,x')+ 

J dx"dp"^^ [u(p', x',p", x")p" 2 t(p",p, x"; z) - u(p', x',po,x")p 2 t(p ,p, x"; z)] 
Pq — p" A 

-i2/jTTp / dx"u(p',x',p ,x") t(p ,p,x";z) , (18) 



where the unknown t(po,p, x"; z) satisfies 
t(p ,p, x';z) = v(p ,p,x')- 



J dx"dp" 2 2/X //2 [u(p , x',p", x")p" 2 t(p" ,p, x"; z) - u(p , x',p , x")p 2 t{p ,p, x"; z)] 
Po P 

-i2fip n / dx"u(p ,x' ,p ,x") t(p ,p,x" ,z). (19) 



These equations define a non-singular set of coupled integral equations. The integrals are approximated by sums 
over quadrature points and weights. The quadrature points are chosen to not include the point po, which is treated 
separately. The resulting linear system is given by 

t(p'i,P, x\;z) = v{pi,p, x\) + 

dx'^dp" 2 2M „ 2 [u(p'i, x' l7 p", x'Dpftijp'! ,p, x'l; z) - u(^,^,p , x'l)plt(p ,p, x'l; z)\ 
kj P ° P J 

~i2irfj,p ^ dx'lu(p', x',p , x'l)t(p ,p, x'l; z) (20) 

k 

and 

t(po,P,x'^ z ) = v (P0,P' X 'l)+ 

X! dx 'k d Pj 2 2M m [ u (Po, x'up", x'l)p" 2 t(p",p, x'l;z) - u(p , x'iiPo, x'l, )plt(Po,P, x'k> z)] 

kj p " Pj 

-i2tthpo ^ dx'lu{p , x' h p , x'l)t(p ,p, x'l; z). (21) 
fe 

This linear system gives approximate solutions for t(p[,p, x\; z) at the quadrature points and at the point where 
Pi =Po- 

Given the approximate solutions at the quadrature points, the results for any values of p',p, x' can be computed 
by inserting the solutions at the quadrature points back in the integral equation (l6| . 

t(p',p,x';z) = v(p',p,x')+ 

y2 dx'ldp" 2 2/i m [u(p\ x',p", x'l)p" 2 t(p",p, x'l; z) - u(p', x',p , x'l)plt(p ,p, x'l; z)] 
jk >'■' 

-i2irfipo ^2 dx'lu{p', x',p , x'l)t(p ,p, x'l; z). (22) 
fc 

These equations give the off-shell transition matrix elements; half-shell matrix elements are obtained by setting 
p 2 = Pq, while on-shell matrix elements have p 2 = p' 2 = p\. 

When calculating NN observables, it is sufficient to solve for the half-shell transition matrix. This solution can be 
used as input 0, E3 , using the first resolvent equation [l8[ , to find the off-shell transition matrix. For this reason we 
only discuss the solution of the half-shell equation, p 2 = p\. 

For realistic models the interactions and transition operators in the above equations are matrices in the spin and 
isospin degrees of freedom. To treat this we expand t(p'j,p, x' k ; z) as a linear combination of known spin operators 
and unknown coefficients functions. A similar expansion is used to represent the interaction except the coefficient 
functions are known for the interactions. The resulting equations arc for the unknown coefficient functions. The 
explicit form of the algebraic equations for the transition matrix elements are given in Appendix [X] 
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IV. CALCULATION OF THE WOLFENSTEIN AMPLITUDES 
A. Formal Considerations 

As noted in Ref. [Hj], one of the major challenges in the calculation of NN observables is the calculation of the 
on-shell amplitudes. In order to avoid taking on-shell limits of off-shell amplitudes expressed as linear combinations 
of six independent spin-operators, we directly calculate the on-shell transition amplitudes by expanding them using 
the five linearly independent Wolfcnstein operators given by Eq. ([?])■ The spin-spin operator in the NN-potential is 
on-shell linearly dependent on the Wolfenstein operators via Eq. (fl4j). since here R = 0, 



(oyK) (<t 2 -K) (a, ■ Q) (<x 2 ■ Q) (er 1 ■ N); (<t 2 ■ N) 

<T1-<T2 = J~ 2 + ^ + ^ , (23) 

where we used the unnormalized vectors, K, Q, and N. The normalization factors are 

K 2 = 2p 2 {l-x) 1 Q 2 = 2p 2 Q (l+x), N 2 = p 4 (l-x 2 ), (24) 

with x := p ■ p' and the magnitude of the on-shell momentum po- For numerical computations it is essential to avoid 
zeros in a denominator. Thus we re-express the operator (<ri • N)(er 2 • N) as 

(<ri • N) (<r 2 ■ N) = 

pi(l x 2 ) a 1 ^ 2 - Pl{l 2 X) (a 1 ■ Q)(<7 2 • Q) - ^ (1 2 +l) ( ffl • K)(er 2 • K), (25) 

which is well behaved in the limits x — > +/ — 1 and po — > 0. This justifies our choice to use Eq. (|25|) to replace the 
operator (<Tx • N)(er 2 • N) in the Wolfenstein basis with the spin-spin operator W 2 = <J\ ■ cr^- 

For the sixth operator we considered two choices: One of them is the time-odd operator W 6 = ((<Ti • K) (er 2 • Q) + 
(<7i • Q)(er 2 ■ K)). The coefficient function in the NN potential for this operator vanishes on-shell. 

A second choice is the Wolfenstein operator (<Xi • N) (<r 2 • N) This operator is related to the time-odd operator, 
W 6 , by re-expressing Eq. (|14p as 

(<T 1 .N)(lT2-N) = 

^Vi • <ra) + (<ri • K)(o- 2 • K) + ((<n • K)(<r 2 ■ Q) + (a 1 ■ Q)(<x 2 • K)) 
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-2 



^- I • Q)(<t 2 • Q). (26) 

If we choose the sixth operator to be (cri ■ N) (<r 2 • N), the on- and off-shell potentials coefficient functions are related 
by 

yonshell Tyoffshcll 

pronshell = yoffshell _ p 4 ^ _ ^proffshell 
yonshcll _ yoffsholl _ Ppi^ + X ) yoffskell 

T/onsholl _ TAoffsholl PpOj ~ X ) T /offshell 

yonshell ^offshcll (27) 

where are the expansion functions for the operators i = 1,2,3,4,5. These expressions show a rapid variation of 
these coefficients near the on-shell point, making the direct computation of the on-shell values in Eq. (|22|) essentially 
impossible. 

However, when the time-odd operator, W 6 , is chosen as 6th operator, the relationship between the on-shell and 
off-shell coefficient functions is 

^Aonshcll proffshcll 

^ronshcll ^offshcll 

-^onshcll -^offshcll 

^Aonshcll proffshcll 

yonshell ^offshcll (28) 

which does not exhibit rapid variations near the on-shell point. Thus, calculations with the time-odd operator as part 
of the basis have the advantage that the coefficient functions are on- and off-shell the same. Test calculations using 



8 



(<Ti • N) (er 2 • N) as 6th operator led to numerical instabilities in solving the Lippmann-Schwingcr equation, as already 
pointed out in Ref. [13| . Thus, our calculations use the time-odd operator, W 6 , as the 6th operator in the basis. 
The scattering amplitude matrix can be expressed in terms of the on-shell transition matrix elements as 

5 

M(p',p) = -4m 2 fd{po,Po,x,z) = -4:Tr 2 ^^2ti{po,p ,x,z)W i (po,po,x), (29) 

i=l 

where we do not display the spin-isospin parameters. 

Accounting for all of the invariances, M can be expressed in the basis of operators of Eq. © in terms of the 
Wolfcnstcin parameters {a, c, m, g, h}. These complex coefficients encode all of the information contained in the on- 
shcll transition matrix elements. By writing the Wolfenstein operators in term of the operators of Eq. (| 1 1[) . the matrix 
M can be expressed in terms of the Wolfenstein parameters as 

M = aW 1 + mW 2 + (g - h - m)W 3 + (g + h - m)W 4 - icW 5 . (30) 

For identical nuclcons this expression has to be anti-symmetric under the exchange of the nucleons, including the 
isospin factors. This can be achieved by anti-symmetrizing the initial two-nucleon state. In our operator formalism 
a particle interchange is represented by a spin exchange and the reversal of the initial momentum, p — > p. The 
spin-exchange operator is given by 

Pl2 = 1 -(l + a 1 -a 2 ) = 1 -(l + W 2 ). (31) 

If we change the sign of the initial momentum p, the vectors K, Q, and N become 

K -> Q 
Q -> K 

N -> -N, (32) 

which implies that the spin-basis elements transform as 

W 1 -> W 1 
W 2 -> W 2 
W 3 -> W 4 

w A -> w 3 

W 5 -> -W 5 . (33) 
The coefficients of the exchange matrix elements become 

w 1 -> w 1 p 12 = + W 2 ) = - (W 1 + W 2 ) 

W 2 -> W 2 P 12 = -jW 2 {I + W 2 ) = TjiZW 1 - W 2 ) 

1 o 2 
W 3 ->■ W 4 P 12 = ^W 4 {I + W 2 ) = w 4 + -^-(w 1 - W 2 ) 

1 K 2 
W 4 -> W 3 P 12 = -W 3 {I + W 2 ) = W 3 + — {W 1 - W 2 ) 

W 5 ->• -W 5 P 12 = -^W 5 {I + W 2 ) = -W 5 . (34) 

The expansion coefficients of the W l become 

ti(po,po,x) ->• U(p ,po, -x). (35) 

Putting all of this together using t = 1 to denote the symmetric iso-triplet subspace, t = to denote the antisymmetric 
iso-singlet subspace, mj(p',p) := —An 2 [i,ti(po,po,x) and rrii(p',— p) = — Air 2 [iti(po,po, — x) we get the following 
expansion for scattering amplitude matrix including exchange contributions 

M = W 1 (mi(p',p) + i(-l)*mi(p', -p) + |(-l)*m a (p', -p)+ 

+ i(-l)*m 3 (p', -p) + ^(-lYm^p', -p) 
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+W 2 (m 2 (p',p) + I(-i)* mi (p', -p) - I(-l)* m2 (p', -P) 
-I(-l)*m 3 (p', -p) - I(-i)* m4 (p', _ P )^ 

+W 3 (to 3 (p',p) + (-)*TO4(p',-p)) 
+W A (m 4 (p / ! p) + (-) t m 3 (p / ,-p)) 

+W 5 (m 5 (p', p) - (-) t m 5 (p' J -p)) . (36) 

Comparing this with Eq. (|30|) gives the following expressions for the Wolfenstein parameters 

1 3 
a = mi(p', p) + -(-l)* mi (p', -p) + -(-l)'m 2 (p', -p) 

+ i(-l)*m 3 (p', -p) + i(-l) t m 4 (p', -p) 

m = m 2 (p', p) + (-l)*im 1 (p / , -p) - i(-l) t m 2 (p' ! -p) 

-I(-l)*m 3 (p', -p) - i(-l)*m 4 (p', -p) 
c = i (m 5 (p', p) - (-)*m 5 (p', -p)) 
g + h = m 2 (p', p) + m 4 (p', p) + --(-l)*m 2 (p', -p) 

+^(-) t m 3 ( P / , -P) - i(-l)*m 4 ( P ', - P ) + ^(-l/milp', "P) 

g-h = m 2 (p',p) + m 3 (p',p) - -(-l) t m 2 (p' ) -p) 

+5(-) t m 4 (p' J -p) - i(-l)*m 3 (p', -P) + ^(-l)'™^', -P). (37) 



B. Wolfenstein Amplitudes for the AV18 Potential 



In this subsection we calculate the Wolfenstein amplitudes for neutron-proton (np) and proton-proton (pp) scat- 
tering calculated with our 3D formulation at 100 MeV and at 350 McV laboratory energy and compare them with 
calculations based on summing partial wave amplitudes up to a given j max . We chose 100 MeV as the low energy, 
since here already the sum over a small number of partial waves should suffice to achieve agreement with the 3D 
calculation. This is indeed the case as illustrated in Fig. [1] for the np Wolfenstein amplitudes. The solid line repre- 
sents the 3D calculations, which perfectly agrees with a partial-wave calculation in which j max = 6. The figure also 
indicates that at this energy even j max = 4 already converges to the 3D result. We also compare our calculations to an 
extraction of the np Wolfenstein amplitudes from the current solution of the GW-DAC Data Analysis Center (2lT[23T ] . 

In Fig. [2] the Wolfenstein amplitudes for pp scattering are shown for 100 MeV laboratory scattering energy, and 
again the 3D calculation is compared to a partial wave sum up to j ma x = 6. Here pp means that we use the strong 
pp interaction as given by the AV18 interaction, however do not treat the long-range Coulomb interaction. The 3D 
and partial-wave summed calculation agree rather well with each other. Only the Real c amplitude displays a small 
deviation. However, when considering the scale of the y-axis, this deviation should not be considered significant. 
In contrast to the np case, the real part of the amplitudes g and h indicate that for very small angles a j max = 4 
calculation is not yet fully converged. 

Figs. [3] and 0] show the Wolfenstein amplitudes for np and pp scattering at 350 MeV laboratory kinetic energy, the 
highest energy for which the AV18 interaction is fitted to the NN data base. In both figures the solid line represents 
the 3D calculation, whereas the dashed line represents a partial- wave calculation summed up to j max = 9. In Fig. [3] 
the partial-wave sums up to j ma x = 2 and jmax = 4 are also shown. In contrast to the lower energy of 100 MeV in 
Fig- HI a t this energy the partial- wave sum up to jmax = 4 is clearly not converged. In fact, the backward angles of the 
real parts of a and h show a particularly bad conv erg ence. In addition, we indicate the extraction of the Wolfenstein 
amplitudes from the GW-DAC data analysis [2l|, HH by the filled diamonds. For the pp Wolfenstein amplitudes in 
Fig. [4] we also display a partial wave calculation summed up to j max = 6 to indicate that even this partial wave sum 
is not yet converged to the full 3D calculation. A careful inspection of the real parts of the forward and backward 
angles of a, <?, and h reveals that even the partial wave sum of j max = 9 is not quite converged yet. In contrast, the 
imaginary parts of all amplitudes show already convergence at roughly j max = 4. 
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C. Numerical Details 



The LS integral equation that needs to be solved when performing a calculation without employing a partial- 
wave decomposition is, in the form of Eqs. (|18|) and (|19|) . a two-dimensional integral equation in two variables (the 
magnitude of a momentum p' and an angle x'). Due to the structure of the equation, an angle integration over the 
azimuthal angle </> can be carried out separately as given in Eq. (|17|) . and is thus not part of the integral equation. 
For the integration we use standard Gauss-Legendre integration with a tangent map of the points of the momentum 
integration. Discretized the set of coupled equations, Eqs. (fT8| and (fT9| . turn into a finite set of linear algebraic 
equations, which is solved by standard methods. 

However, we need to test the numerical convergence for each integration. In the following we demonstrate the 
numerical convergence by concentrating on the real part of the Wolfenstein a amplitude and consider its values for 
different sets of grid points of p, x = cos#, and <fi. For the choice of n p-points, m x-points and I (^-points let us define 
the value of the real part of the Wolfenstein amplitude A as a(n, m, I). Then we define as logarithmic error in cf>, 



Log-Error (7) = log 



(a(60, 30, 40) - a(60, 30, 1)) * 100 



a(60,30,40) 



(38) 



This Log- Error is plotted in Fig. [5jc) for values I £ (5,35), and we see that the <fi integration converges very quickly. 
Already 20 points are sufficient to achieve an error in the order of 10 _4 %. For most of our calculations we choose 
30 4> points. For the x-points wc define 



Log — Error^, (m) = log 



(o(60, 45, 20) - a(60, m, 20)) * 100 



a(60,45,20) 



(39) 



This Log-Error is plotted in Fig. [SJb) for values m 6 (10,40). The x-integration converges much slower with respect 
to the number of mesh-points, and for having an error in the order of 10 _3 to 10 _4 % we need at least 40 integration 
points for a converged calculation. 

For the momentum-integration points wc define 



Log — Error p (n) = log 



(a(75, 30, 20) - a(n, 30, 20)) * 100 



a(75,30,20) 



(40) 



This Log-Error is plotted in Fig.(5ja) for values n € (20, 70). Here we see that the error in the momentum integration 
decreases slowest as a function of the momentum-integration points, and in order to have an error of about 10 _2 % and 
we choose 60 p-points for our calculations. Based on the above considerations, all of our calculations use (n, I, m) = 
(30, 40, 60) mesh points to achieve an overall numerical errors of 0.01%, which is well below the accuracy with which 
NN observables are given. 



V. OBSERVABLES 

One of the main purposes of this paper it to test the momentum-space potential AV18 potential given in 0. There 
it was demonstrated that the potential reproduced the deuteron binding energy and the wave functions. However, 
these tests were only sensitive to the np part of the potential. Now we have the opportunity to test the np as well as 
the pp pieces of the AV18 potential given in Ref. In comparing the results of the scattering calculations to data 
we found that the grid for the Fourier transform of the charge symmetry breaking term, which is the 18-th operator 
in Ref. [l| needed to be extended beyond 100 fm _1 , the standard value for all other terms of the AV18 potential. 
Specifically, for the 18th operator the Chebyshev expansion of the Fourier transform is extended to 250 fin -1 . 

The scattering amplitude matrix is related to the on-shell solution of the Lippmann-Schwinger equation as indicated 
in Eq. (|29j) . Having obtained Af(p',p), the spin averaged unpolarizcd differential cross section is given by summing 
over the final spins and averaging over the initial spins as 

% = - 4 Tr(MM% (41) 



whereas general spin observables have the form 



TrjMAWE) 
{U) Tr(MAft) ' 
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Here the operator A is associated with a measurement of the initial spins and the operator B is associated with a 
measurement of the final spin state. The traces can be expressed as homogeneous polynomials of degree two in terms 
of the Wolfcnstcin parameters. 

To test the potential and the computational method we compare calculated np observable as several energies and 
compare them to the observables in the GW-DAC analysis [U [23[ as well as calculations based on a partial-wave 
expansion (20j . For the pp observables we can only make a comparison with the partial- wave based calculations, since 
we do not include the Coulomb interaction. 

Fig. |6] shows the differential cross sections for pp and np scattering at 100, 300, and 500 McV projectile laboratory 
kinetic energies. The solid lines are the result of solving the equations of Sec. Ill with the momentum-space AV18 
potential of ref. @. Both calculations are compared to partial- wave calculations using a partial- wave sum up to 
jmax = 9. The black diamonds arc the np cross sections taken from the GW-DAC phase shift analysis [Il],[23j]. For 
the higher projectile laboratory kinetic energies there is a clear deviation at the very forward and backward angles 
between the 3D calculation and the partial-wave sum, indicating that even the sum to j m ax = 9 is not yet converged. 

Figs. [7] and [8] show selected spin observables for pp and np scattering at 100 and 300 MeV. The following observables 
are shown, 

Tr(MMttr) , . 

P = — (43) 

Tr(MAft) 1 ' 

= TrCM^-NQMy-N)) 

Tr(MMt) 1 ' 

= TV(AJ((r.(NxQpt( g .Q)) 

Tr(MMt) ■ 1 ' 

Here cr stands for cither <j\ or <x 2 depending on label of the particles. The solid lines represent the 3D calculation with 
the the momentum-space AV18 potential of ref. @, whereas the dashed lines represent a partial- wave calculation 
with partial waves summed up to jmax = 9. The diamonds represent the observables taken from the GW-DAC 
analysis [Il],[23|. At 100 MeV projectile kinetic energy there is perfect agreement between the 3D and the partial- 
wave calculations, as is expected. At 300 MeV there are small deviations in the backward angles of D and R' , however 
not as large as for the differential cross section. 



VI. SUMMARY AND CONCLUSIONS 



We formulated and numerically illustrated an approach to treat NN scattering directly with momentum vectors 
together with spin-momentum operators multiplied with scalar function that are only functions of those momentum 
vectors. Our formulation differs from the one suggested previously in Ref. (l3l | in the choice of the operator expansion. 
The basic foundation for an operator expansion rests on the fact that the most general form of a nucleon-nucleon 
interaction can be represented as a linear combination of six linearly independent spin- momentum operators with scalar 
coefficient functions. The Wolfenstein decomposition of the NN scattering amplitude into five linearly independent 
operators is dictated by physical symmetries [14j , while a sixth operator with those symmetry properties does not exist 
on-shell. The choice of this sixth operator is not unique, and different NN potential employ different choices (Til |24|. 

In this work we solve the Lippmann Schwingcr equation by using different numbers of spin-momentum basis elements 
to represent the on- and half-off shell transition matrix elements. We show that a numerically stable discretization 
of the Lippmann-Schwinger equation is obtained for a certain pair of on and off-shell spin bases. The stable choice 
includes a sixth time-odd operator in the off-shell basis. The choice of this sixth operator is the main difference 
between our formulation and calculation compared to the ones in Ref. [l3|. Our calculations are able to directly 
calculate the on-shell elements of the transition amplitude, which was not easily possible with the choice of Ref. (l3j |. 
We should point out that even if a time-odd operator appears in the expansion of the potential, as it does in the 
AV18 potential, the potential itself remains invariant with respect to time reversal, since the expansion coefficients of 
this operator also contain time-odd functions. The scalar coefficient function time-odd operator we consider vanishes 
on shell. 

With those six operators the Lippmann-Schwinger equation becomes a six channel integral equation in two vari- 
ables. The solution is represented by six complex amplitudes in four variables, which on-shell reduce to five complex 
amplitudes in two variables, which are simple linear combinations of the five Wolfcnstcin amplitudes. We calculated 
the Wolfcnstcin amplitudes for np and pp scattering using the momentum-space representation of the AVI 8 potential 
given in Ref. Q for several laboratory projectile kinetic energies and compared our calculations to the results of 
traditional partial-wave calculations summed up to a given j max and verified for low energies that our 3D calculation 
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matches the partial-wave sum exactly, while at even moderately high energies (like e.g. 300 MeV) a partial-wave sum 
with jmax < 9 is not sufficient for a converged partial-wave calculation. 

Appendix A: Transition Matrix Elements 

In order to keep the notation as simple as possible we use shorthand notations for the terms in Eqs. (|20j|21j) : 



hi 
toi 
vu 

Vol 

jk 

Ok 



UU 

uu 

Uoi 



9j 



t(Pi,Po,x'i';p ) 

v(p'i,P,x't) 
v{po,p,x[) 

u(po,x' h p",x'l) 

u(po,x'i,p ,x'k) 
2/i 

Pi - pf + «0+ ' 



(Al) 



In this notation Eqs. (|20II21[) become the linear system, 

tu = vu + E dx'{.dp" ' g ] [uu-jk{p'j) 2 t k - uu-flkPotok] 
kj 

hi = vol + E dx'ldp"gj[u a i tjk (p") 2 t. ]k - u i-okPlt k] 
kj 



il-Kiipo E dx'luu-oktok 

k 

- i2ir/j,p E dx k u i-oktok- 



(A2) 



The interaction and transition matrix elements are represented by linear combinations of the operators, W 3 . Since 
the W 3 operators depend on the momenta, there are different expansions at each quadrature point. We use the 
notation W 3 for the operators used to expand the on-shell matrix elements and W 3 for the operators used to expand 
the off-shell matrix elements, even though they are identical for j = 1, • • • , 5. Lower indices are used to denote the 
momentum variables that appear in the operator expressions. 



tu ~ 



Substituting these expansions into Eqs 

6 6 



Uol;jk '■ = 
U0l;0k '• = 

2| leads to 

E = E v « w « + 



m— 1 
5 

E « m 

71—1 

E v u w ™ 

71—1 

5 

Em ti jTfi 
v 0l W l 
m—1 
6 

uu-m ■■ E <m w w 

m—1 

6 

u U0k vv ilk 

m—1 
6 

E 



toi 



v,i 



vol 



uu-ok 



m—1 
5 

E 



0l0k vv ilk 



Hk vv lk 



(A3) 



m—1 
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and 



J2 dx 'k d p. 



"9i 



6 5 



.m,n— 1 

6 5 



E dx 'k E E «™o aw^wt 



m—1 n—1 



(A4) 



E w m = E v oi^r+ 



Y dx 'l d P. 



m—1 
6 



'ho 



kj 



6 5 



E <ijk(P?)t%W% k W? k - £ £<0) fc Po«W 



m,n— 1 
5 



m—1 n—1 



(A5) 



In order to obtain a linear system for the coefficient functions, t™ and £q™, we first multiply the first equation from 
the left by W™ and the second one by W™ and then take traces over the spins. The resulting linear system is given 
by 



and 



m—1 



Y dx k d P'j9j 



kj 



6 5 



E u iijk(Pj 2 )t?k^{WF'W$ k W? k ) - Y Tt^okPotWWfW^WZ) 



m,n— 1 

6 5 



m—1 n—1 



(A6) 



Y tmwr'wn = Y <Mwr'wry 



Y dx 'k' d p'. 



fa 



kj 



6 5 



E u °ijk(p" 2 WkT*(frr'wi? k wp k ) - YY. u oio k p 2 otz k Tr(wr'w%wn 



.m,n— 1 
5 



(A7) 



The input to Eqs. (|A6|) through (1A7[) involves a large number of traces over the spin operators. To treat these we use 
the symbolic tool developed in PJJ] to evaluate these traces. 
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Wolfenstein Amplitudes np E lab = 100 MeV 
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cm. L ° cm. L ° ' 

FIG. 1. The Wolfenstein amplitudes for neutron-proton scattering at 100 MeV laboratory kinetic energy based on the AV18 
potential. The solid (red) line represents the 3D calculation, whereas the dashed (blue) curve is obtained from a partial-wave 
calculation summing partial waves up to j = 6. The partial- wave sums up to j — 2 and j = 4 are shown as double-dash-dotted 
(turquoise) and dotted (green) lines. The data points are the amplitudes extracted from the GW-DAC current analysis [2ll.[23|. 
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FIG. 2. The Wolfenstein amplitudes for proton-proton scattering a 100 MeV laboratory kinetic energy based on the AV18 
potential. The meaning of the curves is the same as in Fig. [T] 
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Wolfenstein Amplitudes np E lab = 350 MeV 
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FIG. 3. The Wolfenstein amplitudes for neutron-proton scattering at 350 MeV laboratory kinetic energy based on the AV18 
potential. The solid (red) line represents the 3D calculation, whereas the dashed (blue) curve is obtained from a partial-wave 
calculation summing partial waves up to j = 9. The partial-wave sums up to j = 2 and j = 4 are shown as double-dash-dotted 
(turquoise) and dotted (green) lines. The data points are the amplitudes extracted from the GW-DAC current analysis [2ll.[23|. 
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FIG. 4. The Wolfenstein amplitudes for proton-proton scattering at 350 MeV laboratory kinetic energy based on the 
potential. The meaning of the curves is the same as in Fig. [3] 
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np Scattering at E =100 MeV 
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FIG. 5. The logarithmic error for the momentum grid points (a), the polar angle (cos# = x) points (b) and the azimuthal 
angle (cj>) grid points (c) as function of number of points. Precise definitions and explanations are given in Sec. IIV CI 
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FIG. 6. The difTerential cross section for pp (left) and np (right) scattering as function of the cm. scattering angle for 100, 300, 
and 500 MeV projectile laboratory kinetic energy. The solid (red) line represent the 3D calculation with the AV18 potential, 
whereas the dashed line represent a partial-wave based calculation summed up to j ma x = 9. The diamonds represent the np 
data from the GW-DAC analysis [HI, HI). 
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FIG. 7. The spin observables P, D and R' for pp (left) and np (right) scattering as function of the cm. scattering angle 
for 100 MeV projectile laboratory kinetic energy. The solid (red) line represent the 3D calculation with the AV18 potential, 
whereas the dashed line represent a partial-wave based calculation summed up to jmax = 9. The diamonds represent the np 
data from the GW-DAC analysis [MUl]. 




FIG. 8. The spin observables P, D and R' for pp (left) and np (right) scattering as function of the cm. scattering angle for 
300 MeV projectile laboratory kinetic energy. The meaning of the curves is the same as in Fig. [7] 



